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Abstract 

We consider an infinite graphene geometry where bonds and sites have been removed selectively 
to map it onto an effective Sierpinski gasket comprising of hexagons. We show that such a struc- 
ture is capable of sustaining an infinite number of extended single particle states inspite of the 
absence of any translational order. When each basic hexagonal plaquette in the Sierpinski geom- 
etry is threaded by a magnetic flux, the spectrum exhibits bands of extended eigenstates. The 
bands persist for any arbitrary value of the flux but disappear again as the flux becomes equal 
to half the fundamental flux quantum. The localization - de-localization issues are discussed 
thoroughly along with the computation of two terminal magneto-transport of finite versions of 
the lattice. The numerical results corroborate our analytical findings. 
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1. Introduction 



Graphene systems have been at the center stage of research in low dimensional systems, in 
particular, nano-structures, in the last decade. Both graphene systems in 2D and graphene 
nano ribbons have been proposed as potential candidates for various applications in nanotech- 
nology [I] - [7]. A wide variety of problems have been investigated not only from the angle 
of applications, but also from the standpoint of the novel physics offered by such systems. 
Electronic and magnetic properties of triangular graphene rings [8], structural, mechanical and 
electrical properties of defect-patterned graphene nano-meshes [9], molecular transistors con- 
nected to graphene nano-ribbon electrodes [TU] are among some of the recent exciting studies in 
the graphene system. 

In recent years, several interesting experiments and theoretical works have been reported which 
deal with vacancies in a graphene sheet or nanoribbon. Such vacancies are usually created by 
ion bombardment, and in principle, the locations can be controlled almost at will. Pereira and 
Schulz [11] have studied the effects of vacancies on the electronic properties of a graphene sheet 
in the presence of a perpendicular magnetic field. The effect of an extended linear defect on the 
electronic transport properties of a graphene sheet has been investigated by Bahamon, Pereira 
and Schulz [12] . It is shown that such defects profoundly modify the properties of a nano-ribbon 
and even introduce new conductance quantization values. Using a first principles method Jippo, 
Ohfuchi and Kaneta |13j report the transport properties of graphene sheets having two- and 
one-dimensional periodic nano holes. 

Vacancies in graphene lattice are shown to give rise to extra localized states between Landau 
levels [11]. The distance between the vacancies and their pattern of distribution have non-trivial 
effect on the physical properties of a 2D graphene geometry. The modification of the bands, the 
nature of the wave functions in the presence of the vacancies are inspiring in introducing such 
vacancies on purpose following a tailor made design [11], and look for any unusual properties 
that might suggest new applications. 

It is thus intriguing to know what happens to the electronic states of a graphene sheet when 
one generates holes of increasing size, and distributes them over an infinite two dimensional 
graphene sheet following a pre-determined geometry. This is the central motivation behind the 
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present work. We undertake a detailed investigation of the single particle states on an infinite 
graphene sheet where holes with increasing sizes have been distributed following a Sierpinski 
fractal geometry [H] - [16]. 

Such a study serves a dual purpose. Sierpinski gasket (SPG) fractals have been well known 
candidates for studying percolation clusters. These self similar lattices are known to give rise 
to exotic electronic energy spectrum, revealing a fragmented Cantor set character [T3] - |16j . 
that has been shown later to contain an infinite number of extended eigenstates embedded in 
it [T7] - [2D]. Recently, a hexagonal Sierpinski structure has been experimentally synthesized [21] . 
The graphene with a Sierpinski distribution of holes lies very close to the synthesized structure, 
and gives us an opportunity to analyze the character of electronic states in these systems. 
Secondly, the present day lithographic techniques allows us to fabricate networks following any 
desired geometry. The potential of a graphene-fractal in the field of nano devices is thus worth 
investigating. 

We work within a single band tight binding formalism. Several works within the tight binding 
formalism and with spinless, non-interacting electrons have been successful in highlighting the 
physical properties of graphene systems with or without defects [11] - [13] . We use the same, 
together with the real space renormalization (RSRG) group decimation scheme that exploits the 
self-similarity of a fractal geometry, to unravel the character of single particle eigenstates. The 
influence of a magnetic flux piercing a selected subset of the graphene hexagons on its electronic 
properties, is studied in details. In addition to this, the two terminal electronic transport across 
finite but arbitrarily large graphene-fractals is worked out. 

Our results are interesting. We find that in the absence of any magnetic field, the graphene- 
fractal gives rise to an infinite number of extended eigenstates which coexist with the fragmented 
spectrum of localized states. The energies of the extended states can be precisely determined 
from the RSRG recursion relations. As soon as a magnetic field is switched 'on' in an elemen- 
tary hexagon of the graphene-fractal, the density of states exhibits continuous distribution of 
extended states in the spectrum. These bands of extended states persist for all non- zero val- 
ues of the magnetic fiux except for ^ = ^o/2, where <&o = hc/e is the fiux quantum. At 
<^ = <&o/2 the spectrum swings back to the flux free shape with sharply localized eigenstates 
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dominating the entire spectrum. Incidentally, appearance of such bands of extended states has 
been suggested earlier, based on extensive numerical studies by Chakrabarti ^9] and Schwalm 
& Moritz [22] in the context of a 3-simplex fractal and a modified rectangle lattice respectively. 
These are self similar lattices without any translational invariance. The existence of flux driven 
bands of extended states is thus curious. With the results presented in this communication we 
have a strong case where we can conjecture that bands of extended states might be a generic 
feature of deterministic fractals with holes. 

In what follows, we describe the results of our calculation. Section 2 describes the model and 
the method. In section 3 we extract and analyze the extended eigenstates, discuss the density 
of states with and without the magnetic flux threading the basic plaquettes and, present the 
detailed results of the magneto-transport calculations for finite graphene-fractals . In section 4 
we draw the conclusions. 

2. The model and the method 

Let us refer to Fig. [TJ A selective removal of sites from the two dimensional graphene sheet will 
lead to a set of extended defects (voids) which distribute themselves on an effective Sierpinski 
fractal network |14j . With both the A (black) and B (red) sites now present in the lattice, the 
two dimensional geometry in fact resembles a Zachariascn fractal glass [23]. The fluctuating 
environment around each site generates unusual eigenstates and strange transport properties as 
studied earlier [T7] - [20]. We shall be concerned with the effect of the parent graphene geometry 
on the electronic states of such a system with multiple holes. To test the effect of a magnetic 
flux on the spectral properties within a minimal model, we include a uniform magnetic flux $ 
through each of the hexagons in Fig. [TJ The system is described by the standard tight binding 
Hamiltonian, 

H = Y.'^\^)(^\ + E [ U,e^m{j\ + i,.e-^-U)(.| ] (1) 

where, is the on-site potential at the i-th atomic site and tij is the nearest neighbor hopping 
integral. Oij = 27r$/6<&o is the Peierl's phase which is included in the hopping along each arm 
of a hexagon. $o = hc/e is the fundamental flux quantum. With the selective inclusion of the 
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magnetic flux, the graphene- fractal now describes a system where the time reversal symmetry 
is broken only on a subset of the bonds. Consequent changes in the energy spectrum and 
transport characteristics thus need to be examined in details. It may be mentioned that, the 
triangular symmetry underlying a graphene topology has been exploited previously by Vargas 
and Naumis [53] to explain an increase in localization in a doped graphene. The present system 
is thus worth examining. 




Figure 1: Schematic diagram of a portion of the graphene-fractal scooped out of an infinite 
graphene sheet. Each of the surviving hexagonal plaquettes is threaded by a uniform 
magnetic flux The underlying dotted line shows the Seirpinski gasket fractal geometry 
in which the basic graphene-fractal geometry can be mapped onto by decimating the 'B' 
type (red) atomic sites. 

For our purpose we shall consider constant values of ej = eo and tij = to throughout the 
calculation. That is, the A- and the B-sites are not distinguished energetically. In the context 
of the graphene sheet they are all carbon atoms. 

The Schrodinger equation for the original graphene network is equivalently cast into a set of 
difference equations, viz, 

{E - ei) i^i = Y^ Uje"''^ V,- (2) 
j 

which is utilized in decimating out a suitable subset of atoms to reduce the infinite graphene- 
fractal to a standard triangular SPG. In this triangular SPG, there is a uniform value of the 
effective on-site potential e, which, of course, is a function of energy now. The nearest neighbor 
hopping integrals assume two different values depending on the phase associated with them. 
The distribution is illustrated in Fig. EJa). The bonds of the inner triangles are associated with 
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the hopping integrals Tf (forward) and t;, (backward) , while those along the outer triangle are 
designated hy tf and tj, respectively. The values of these parameters are given by, 
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Figure 2: (a) Schematic diagram of the effective Seirpinski gasket fractal which is obtained 
by decimating the the 'B' type (red) atomic sites from the original graphene-fractal network. 
The arrows show the direction for forward hopping from one atomic site to another atomic 
site - for the inner triangles it is tj and for outer triangles it is tf. (b) The renormalized 
version of (a). 



The effective SPG network (Fig. [2](a)) is easily renormalized (Fig. [2jb)) and the renormalized 
values of the parameters are, 

e' = e + 2{utf + u*tb) 

Tf = Vtf + W*tb, = Wtf + V*tb 

t'f = vtf + w*tb, t'i, = wtf + v*tb (4) 
6 



where, u = {p + rp*)/V^ v = [q + rs*)/V, w = [s + rq*) /V 

u* = {p* + r*p)/V, V* = {q* + r*s)/V, w* = (s* + r*q)/V 
T> = 1 — rr* 



Here, „ = , ^ (^Hi^ll^A 

p. ^ iE-e)tf ^ ^ {E-e)U + Tftf 

So ' So 
{E - e)Tf + inf ^ ^ ntf 

So ' So 

^* ^ {E-e)n + {TfY ^ ^.^rjh 

So ' So 

So = {E-ef-TfU 



The above set of recursion relations will now be analyzed to see the effects of a fractal distri- 
bution of holes in a graphene sheet with magnetic flux trapped in a selected set of the hexagons. 



3. Results and discussion 

3.1 Local Density of States (LDOS) at the bulk atomic sites 

We have calculated the LDOS at the 'A' type (black) bulk atomic sites for an infinite graphene- 
fractal using the Green's function and a standard real space decimation technique [2^. The 
LDOS is given by, 

(5) 



Poo{E) = lini 



■-lmGoo{E + ir]) 

IT 



where, Goo{E) is the local Green's function at the 'A' type (black) bulk atomic sites. To obtain 
Goo 5 the recursion relations Eq. Q are iterated with a small imaginary part added to the energy 
E until the magnitude of the nearest neighbor hopping goes to zero (or, equivalently, becomes 
less than a small pre-assigned quantity). The on-site potential flows to a fixed point value e*, 
and in this limit Goo = {E + irj — e*)~^. 



3.1.1 The zero flux case 



In the absence of any magnetic flux the LDOS shows a fragmented structure (Fig. [21 top 
panel). Most of these fragmented states are locahzed states, which can be verified by studying 
the flow the hopping integral under successive iterations, keeping the energy at a particular value. 
With any arbitrarily chosen energy E at which the LDOS is non-zero, the hopping integral flows 
to zero under RSRG steps. 

Interestingly, such a graphene-fractal is found to sustain an inflnite number of extended eigen- 
states. These eigenstates coexist with the localized states. This turns out to be a generic feature 
of the underlying triangular Sierpinski geometry [n\ 118] . To extract such extended states we 
need to set the energy of the electron E = e^^^ , where, e^^^ represents the effective on-site poten- 
tial at a bulk atomic site after the ^-th stage of renormalization. A look at the recursion relations 
given in Eq. with $ = will reveal that, setting E = e*-^^ immediately leads to a 'two-cycle' 
fixed point, viz, e^^"*"^^ = e*-^"*"^^ = e^^^ and i^^"*"^) = — = t^^^ beginning at a certain stage i 
of renormalization. E — = gives a polynomial equation of E, the real solutions of which will 
yield the values of energies at which the states are extended. Of course, the solutions will have 
to lie within the spectrum of the graphene-fractal. 

For example, if we set E = e at the basic stage £ = 0, with eo = 0, and to = Ij then E = ±^/2 
will be energies for two extended states sustained by the fractal lattice in the absence of any 
magnetic fiux. If we set E = e^^^ on a one step renormalized lattice and with zero flux, it leads 
to an equation, 

E^ - 8E'^ + 17^2 _ 10 = (6) 

The roots are, E = ±1, ±V2, ±^/5. Out of these ^ = ±1 are 'spurious' roots, not included in 
the spectrum of the graphene-fractal, while the others are there, and correspond to two-cycle 
fixed points of the RSRG recursion relations beginning at the stage i = 1. These are extended 
eigenstates. It is to be noted that the roots evolving from the solution of the equation E = e^^^ 
are also included in the spectrum of E = e^^"^^^. This is again a generic feature of the SPG |18] . 
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Figure 3: Plots for Local Density of states (LDOS) witli the energy of the electron (E) at 
the bulk atomic sites ('A' type) of an infinite graphene-fractal network for different values 
of magnetic flux (a) represents the case for $ = 0, (b) represents the case for <I> = <I?o/4 
and (c) represents the case for $ = ^o/2. We have set eo = and to = 1. 

3.1.2 The non-zero flux cases 

(a) The general character of the spectrum: 

As we set the magnetic flux in each hexagon to a non-zero value, continuous bands are seen 
to be created in the LDOS profile (Fig. [3l middle panel). In this we have shown the LDOS for 
= ^o/4: and it clearly shows the continua in the range 1 < \E\ < 2. In such continua the states 
are 'extended', as has been examined by picking up any energy in this domain arbitrarily, and 
studying the flow of the hopping integral. The hopping integral in general, oscillates chaotically 
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and remain non-zero for an indefinite number of iterations. The LDOS profile shows a continuum 
zone for all non-zero values of < $ < ^o/2- Energy eigenvalues, picked up arbitrarily from 
any such continuous zone, corresponds to an extended eigenstate of the system. The localized 
eigenstates are also there. The spectrum is a mixture of the bands of extended states and 
clusters of localized states. However, it is not apparent whether there is any possibility of a 
metal- insulator transition driven by the magnetic flux. At $ = the fragmented character 
in the LDOS profile is restored again (Fig. [3l bottom panel). 

(b) Fixed points: 

In the presence of a magnetic flux, the extraction of a fixed point of the transformations Eq. 
is a non-trivial issue. Nevertheless, it is possible to work out the case for £ = 1. Setting E = e^^^ 
leads to the equation, 

f{E) = cos 69 (7) 
where, f{E) = [E^ - SE"^ + llE"^ - 8)/2 and 6 = 27r$/6$o- Clearly, for any energy E lying 




-1 1 

Energy (E) 



Figure 4: Plot of energy function f{E) on the LHS of Eq. ([7]) with energy E. The shaded 
portions under the curve shows the values of energy for which the value of f{E) lies within 
±1. 



within the spectrum, whenever the left hand side of Eq. d?]) remains bounded by ±1, one can 
tune the magnetic flux appropriately to satisfy the equation E = e^^\ For this special value of 
the magnetic flux, a two- cycle fixed point behavior is obtained. We depict the variation of the 
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left hand side of Eq. ([7]) in Fig. HI In principle, any energy within the shaded zone (provided it is 
included in the overall energy spectrum) can be made to correspond to an extended eigenfunction 
by a magnetic flux, whose value can be estimated from Eq. ([7]). For example, if we set the Fermi 
level at E = 0.7 (in unit of Iq), then it can be worked out to find that a value of the magnetic 
flux $ = 0.38178$o throws the parameter space into a two-cycle fixed point beginning 

at the i = 1. As soon as one shifts away from this special value of the flux, the hopping integrals 
usually flow to zero as the renormalization progresses. This implies that, once we fix the Fermi 
level at a special value, the entire graphene-fractal can be made to conduct or act as an insulator 
by tuning the magnetic flux piercing the elementary hexagons. 

(c) Flux controlled behavior of the eigenstates: 

Certain points are of interest, and need to be appreciated in the context of flux controlled 
behavior of the eigenstates in a graphene-fractal. To appreciate this discussion we first illustrate 
in Fig. [5] the variation of the LDOS against the magnetic flux at special values of energy 
(Fig. [5]^a)), and for E = ^/b (Fig. [5]^b)). In each diagram it is clear that the selected energy 
can be brought inside or thrown out of the energy spectrum by tuning the magnetic flux. We 
specially focus in and around <1> = 

Let us first consider the case when E = There is clearly a continuous distribution of flux 
values in the neighborhood of $ = <I>o/2, viz, for 0.495<l>o < <I> < 0.505$o, for which the LDOS 
is non-zero. For this range of magnetic flux, with the selected energy, the hopping integral does 
not flow to zero, but becomes small and keeps on oscillating around a value ~ 10~^. This implies 
that the wavefunction really does not decay, but the overlap of amplitudes at distant lattice sites 
becomes very small. This tempts us to categorize such states as critical [T5], or at least they 
have very large localization lengths. 

Choosing E = ^/b leads to a more interesting scenario. We now have a sharply localized 
eigenstate for <I> = <l*o/2. The LDOS is isolated at $ = $o/2, and is highly degenerate. Setting 
E = \/5, and <I> = <l*o/2 makes the hopping matrix element flow to zero quickly. This obser- 
vation speaks in favor of the localized character of the states. Interestingly, in the immediate 
neighborhood of the center of the spectrum ($ = $o/2), the hopping integral displays non- zero 
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Figure 5: Variation of Local Density of States (LDOS) with tlie magnetic flux ^ at the 
bulk atomic sites ('A' type) of an infinite graphene- fractal network, (a) for energy 
and (b) for energy E = ^/b. We have set eo = and to = 1- 



values under successive iterations, indicating either completely extended states, or critical ones 
(in the sense as discussed above). We thus encounter a possibility of a reentrant crossover in 
the nature of the spatial extension of the wavefunction at £' = \/5 by tuning the magnetic flux 
around the half flux quantum. 



3.2 Transmission characteristics of the finite graphene- fractal network 

To get the two terminal end-to-end transmission coefficient for a finite sized graphene-fractal, 
we clamp the system between two semi-infinite ordered leads. The leads, in the tight binding 
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model, are described by a constant on-site potential and a nearest neighbor hopping integral 
tl- We then successively renormalize the system to reduce it into an effective dimer (Fig. [6|) 
consisting of two 'renormalized' atoms, each having an effective on-site potential equal to e and 
with an effective hopping integral t. 
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Figure 6: Schematic diagram of the effective dimer clamped between two semi-infinite 
ordered leads. 



The transmission coefficient across the effective dimer is given by 

^ 4 sin^ ka , . 

T = 5 ^ 8) 

[(Mi2 - M21) + (Mil - M22) cos kaY + [(Mii + M22) sin ka] 

where, 

Mil = ^^TT^ - 7^' = M21 = -M12, M22 = 

tfti tL tf tf 

and 'a' is the lattice constant and is taken to be equal to unity throughout the calculation. 

In Fig. [7] (a), (b), and (c) we plot the two terminal transmission coefficient of a 4-th generation 
graphene-fractal for <1> = 0, for $ = <l>o/4 and ^ = respectively. The figures bring out the 

typical fragmented spectrum for zero magnetic field, which then gets converted into patches of 
continua at $ = ^o/4. Finally, at $ = $o/2, the fragmented character is restored, reflecting a 
poorly conducting system. These results are at par with our discussion about the general band 
structure of the graphene- fractal system. 

We have also examined the AB oscillations in the transmission spectrum at particular values 
of the electron energy E. The oscillations have the typical <I>o periodicity observed in systems 
even with a single loop. However, we do not show this result to save space. 
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Figure 7: Transmission coefficient across a 4-th generation graphene- fractal network for 
different values of magnetic flux ^. (a) represents the case for $ = 0, (b) represents the 
case for $ = ^o/4 and (c) represents the case for $ = ^o/2- We have set eo = cl = and 
to = tL = 1. 



4. Concluding remarks 

In conclusion, we have considered an infinite graphene sheet in which holes have been created 
following a Sierpinski fractal distribution. The electronic spectrum of the system, investigated 
within a tight binding Hamiltonian for spinless, non-interacting electrons reveal a wide variety in 
the nature of single particle states. The magnetic field piercing each basic hexagonal plaquette 
is shown to lead to absolutely continuous parts in the spectrum. The results are reflected in the 
flow of the hopping integrals, as observed within a real space renormalization group formalism. 
A crossover in the behavior of the wavefunction at particular values of the energy can be made 
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to happen by tuning the magnetic flux. The graphene-fractal system is thus a candidate to be 
inspected more carefully as a potential electronic device. 
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